Phase transitions of choline dihydrogen phosphate: A vibrational spectroscopy and periodic DFT study.

Choline dihydrogen phosphate, [Chol][H2PO4], is a proton-conducting ionic plastic crystal exhibiting a complicated sequence of phase transitions. Here, we address the argument in the literature around the thermal properties of [Chol][H2PO4] using Raman and infrared microspectroscopy. The known structure of the low-temperature crystal, which contains the anti-conformer of [Chol]+ and hydrogen-bonded dimers of anions, was used to do periodic density functional theory calculations of the vibrational frequencies. Raman spectra indicate that the solid-solid transition at 20 °C is linked to a conformational change to the gauche [Chol] conformer with a concurrent local rearrangement of the anions. The distinct bands of lattice modes in the low-frequency range of the Raman spectra vanish at the 20 °C transition. Given the ease with which metastable crystals can be produced, Raman mappings demonstrate that a sample of [Chol][H2PO4] at ambient temperature can contain a combination of anti- and gauche conformers. Heating to 120 °C causes continuous changes in the local environment of anions rather than melting as suggested by a recent calorimetric investigation of [Chol][H2PO4]. The monotonic change in vibrational spectra is consistent with earlier observations of a very small entropy of fusion and no abrupt jump in the temperature dependence of ionic conductivity along the phase transitions of [Chol][H2PO4].

DFT and ab initio molecular dynamics simulation study of the infrared spectrum of the protic ionic liquid 2-hydroxyethylammonium formate.

Protic ionic liquids (PILs) typically show a complex band shape in their infrared (IR) spectra in the high-frequency range due to the hydrogen stretching vibrations of functional groups forming rather strong hydrogen bonds (H-bonds). In the low-frequency range, the intermolecular stretching mode of the H-bond leaves a mark in the far-IR spectrum of PILs. In this study, the IR spectrum of the PIL 2-hydroxyethylammonium formate, [HOCH2CH2NH3][HCOO], is investigated in order to identify the different modes that contribute to the high-frequency band shape, i.e. the cation ν(NH), ν(OH), and ν(CH) modes, and the anion ν(CH) mode, as well as the intermolecular mode of the strongest H-bond in the far-IR spectrum. The assignment is validated by quantum chemistry calculations of clusters at the density functional theory (DFT) level for four ionic pairs and by ab initio molecular dynamics (AIMD) simulations of ten ionic pairs. There is good agreement between the vibrational frequencies obtained from DFT and AIMD simulations for both the high- and low-frequency ranges. Based on the calculations, the strong H-bond interaction between the cation -NH3 group and [HCOO]- gives a broad band envelope associated with the ν(NH) mode in the high-frequency range of the IR spectrum on which there are narrower peaks corresponding to the ν(OH) and ν(CH) modes. In the far-IR (FIR) spectrum, the anions' rattling motion gives a broad feature with a maximum at 160 cm-1, while the H-bond's intermolecular NH⋯O stretching mode appears as a peak at 255 cm-1.

Confined ionic liquids films under shear: The importance of the chemical nature of the solid surface.

Ionic liquids have generated interest in applications as lubricants and as additives to conventional lubricants due to their unique physical properties. In these applications, the liquid thin film can be subjected simultaneously to extremely high shear and loads in addition to nanoconfinement effects. Here, we use molecular dynamics simulations with a coarse-grained model to study a nanometric film of an ionic liquid confined between two planar solid surfaces both at equilibrium and at several shear rates. The strength of the interaction between the solid surface and the ions was changed by simulating three different surfaces with enhanced interactions with different ions. The increase in the interaction with either the cation or the anion leads to the formation of a solid-like layer that moves alongside the substrates; however, this layer can exhibit different structures and stability. An increase in the interaction with the high symmetry anion produces a more regular structure that is more resistant to the effects of shear and viscous heating. Two definitions were proposed and used for the calculation of the viscosity: a local definition based on the microscopic characteristics of the liquid and an engineering definition based on the forces measured at the solid surfaces, with the former displaying a correlation with the layered structure induced by the surfaces. Because of the shear thinning behavior of the ionic liquids as well as the temperature rise brought on by viscous heating, both the engineering and the local viscosities decrease as the shear rate increases.

Synthesis and Physicochemical Properties of Acrylate Anion Based Ionic Liquids.

Two polymerizable ionic liquids (or monomeric ionic liquids, mILs) namely 1-butyl-3-methylimidazolium and choline acrylates ([C4mim]A and ChA, respectively) were synthesized using the modified Fukumoto method from corresponding chlorides. The chemical structure of the prepared mILs was confirmed with FTIR and NMR study. Investigation of the thermal properties with DSC demonstrates that both mILs have a Tg temperature of about 180 K and a melting point around 310 K. It was shown that the temperature dependence of FTIR confirm the Tg to be below 200. Both mILs exhibited non-Newtonian shear thinning rheological behavior at shear rates >4 s−1. It was shown that [C4mim]A is able to dissolve bacterial cellulose (BC) leading to a decrease in its degree of polymerization and recrystallisation upon regeneration with water; although in the ChA, the crystalline structure and nanofibrous morphology of BC was preserved. It was demonstrated that the thixotropic and rheological properties of cellulose dispersion in ChA at room temperature makes this system a prospective ink for 3D printing with subsequent UV-curing. The 3D printed filaments based on ChA, containing 2 wt% of BC, and 1% of N,N'-methylenebisacrylamide after radical polymerization induced with 1% 2-hydroxy-2-methylpropiophenone, demonstrated Young's modulus 7.1 ± 1.0 MPa with 1.2 ± 0.1 MPa and 40 ± 5% of strength and ultimate elongation, respectively.

Molecular Dynamics Simulation Study of the Far-Infrared Spectrum of a Deep Eutectic Solvent.

Deep eutectic solvents (DESs) are similar to ionic liquids (IL) in terms of physicochemical properties and technical uses. In ILs, far-infrared (FIR) spectroscopy has been utilized to reveal ionic interactions and even to produce a signature of the strengthening of the cation-anion hydrogen bond. However, for the situation of the DES, where the mixing of a salt and a molecular species makes the interplay between multiple intermolecular interactions even more complex, a full investigation of FIR spectra is still absent. In this work, the FIR spectrum of the DES, often referred to as ethaline, which is a 1:2 mixture of choline chloride and ethylene glycol, is calculated using classical molecular dynamics (MD) simulations and compared to experimental data. To explore the induced dipole effect on the computed FIR spectrum, MD simulations were run with both nonpolarizable and polarizable models. The calculation satisfactorily reproduces the position of the peak at ∼110 cm-1 and the bandwidth seen in the experimental FIR spectrum of ethaline. The MD simulations show that the charge current is the most important contributor to the FIR spectrum, but the cross-correlation between the charge current and dipole reorientation also plays a role in the polarizable model. The dynamics of the chloride-ethylene glycol correlation span a wide frequency range, with a maximum at ∼150 cm-1, but it participates as a direct mechanism only in the charge current-dipole reorientation cross-term. Anion correlations, whose dynamics are regulated via correlation with both ethylene glycol and choline, make the most significant contribution to the charge current mechanism. The MD simulations were also utilized to investigate the effect on the FIR spectrum of adding water to the DES and switching to a 1:1 composition.

Role of density and electrostatic interactions in the viscosity and non-newtonian behavior of ionic liquids - a molecular dynamics study.

Strong ionic interactions, as well as the consequent correlations between cation and anion dynamics, give ionic liquids various physical features that set them apart from ordinary organic solvents. In particular, they result in larger viscosities and larger densities than mixtures of neutral compounds with similar molecular structures. However, both the direct effect of electrostatic interactions and the increase of liquid density contribute to the high viscosity and so far no experimental or computational work enabled a clear quantification of those effects. Also, the effects over the shear thinning behavior, which may have important consequences for application as lubricants, were not considered yet. Here, these questions were tackled by performing non-equilibrium molecular dynamics (NEMD) simulations changing both the strength of ionic interactions and liquid density at several shear rates using a coarse grained model. The relative dielectric constant was adjusted to reproduce viscosity data from all-atoms simulations on both zero shear and high shear conditions. Elimination of ionic interactions results in a reduction of density and zero shear viscosity and also delays the beginning of shear thinning to higher shear rates. Restoring density to the ionic liquid's value only partially reverses the alterations. Correlations of the non-newtonian behavior and changes in the intermolecular structure and contact lifetimes were also explored.

Relating the structure and dynamics of ionic liquids under shear by means of reverse non-equilibrium molecular dynamics simulations.

The effect of the shear rate on the viscosity and the structure of 1-ethyl-3-methylimidazolium based ionic liquids with three different anions (tetrafluoroborate, dicyanamide, and bis(trifluoromethylsulfonyl)imide) was studied by means of reverse non-equilibrium molecular dynamics (RNEMD) simulations using a polarizable force field. The three liquids display a Newtonian plateau followed by a shear thinning regime at shear rates of the order of GHz. Even though the main features of the liquid structure remains under shear, systematic changes were noticed at the GHz rates, with coordination shells becoming more diffuse as noticed by the reduction in the difference between consecutive maxima and minima in the radial distribution function. Interestingly, these structural changes with the shear rate can be precisely fitted using the Carreau equation, which is a well-known expression for the shear rate dependence of the viscosity. The fitting parameters for different distributions can be used to explain qualitatively the shear thinning behavior of these liquids. In the GHz range, the cations and, in a minor extension, some anions, tend to assume preferentially a parallel orientation with the flux, which contributes to the shear thinning behavior and may have consequences for adhesion in applications as lubricants.

Effect of alkyl-group flexibility on the melting point of imidazolium-based ionic liquids.

The low melting point of room temperature ionic liquids is usually explained in terms of the presence of bulky, low-symmetry, and flexible ions, with the first two factors related to the lattice energy while an entropic effect is attributed to the latter. By means of molecular dynamics simulations, the melting points of 1-ethyl-3-methyl-imidazolium hexafluorophosphate and 1-decyl-3-methyl-imidazolium hexafluorophosphate were determined, and the effect of the molecular flexibility over the melting point was explicitly computed by restraining the rotation of dihedral angles in both the solid and the liquid phases. The rotational flexibility over the bond between the ring and the alkyl chain affects the relative ordering of the anions around the cations and results in substantial effects over both the enthalpy and the entropy of melting. For the other dihedral angles of the alkyl group, the contributions are predominantly entropic and an alternating behavior was found. The flexibility of some dihedral angles has negligible effects on the melting point, while others can lead to differences in the melting point as large as 20 K. This alternating behavior is rationalized by the different probabilities of conformation defects in the crystal.

An inelastic neutron scattering, Raman, far-infrared, and molecular dynamics study of the intermolecular dynamics of two ionic liquids.

The intermolecular dynamics in the THz frequency range of the ionic liquids n-butyl-trimethylammonium bis(trifluoromethanesulfonyl)imide, [N1114][NTf2], and methyl-tributylammonium bis(trifluoromethanesulfonyl)imide, [N1444][NTf2], were investigated by a combined usage of inelastic neutron scattering (INS), Raman, and far-infrared (FIR) spectroscopies and the power spectrum calculated by molecular dynamics (MD) simulations. The collective dynamics of the simulated systems is also discussed by the calculation of time correlation functions of charge and mass currents that are projected onto acoustic- and optic-like motions. The INS and Raman measurements have been performed as a function of temperature in the glassy, crystalline, and liquid phases. The excess in the vibrational density of states over the expectation of the Debye theory, the so-called boson peak, is found in the INS and Raman spectra as a peak at ∼2 meV (∼16 cm-1) and also in the direct measurement of heat capacity at very low temperatures (4-20 K). This low-frequency vibration is incorporated into the curve fits of Raman, FIR, and MD data at room temperature. Fits of spectra from these different sources in the range below 100 cm-1 are consistently achieved with three components at ca. 25, 50, and 80 cm-1, but with distinct relative intensities among the different techniques. It is proposed as the collective nature of the lowest-frequency component and the anion-cation intermolecular vibration nature of the highest-frequency component. The MD results indicate that there is no clear distinction between acoustic and optic vibrations in the spectral range investigated in this work for the ionic liquids [N1114][NTf2] and [N1444][NTf2]. The analysis carried out here agrees in part, but not entirely, with other propositions in the literature, mainly from optical Kerr effect (OKE) and FIR spectroscopies, concerning the intermolecular dynamics of ionic liquids.

Pseudo-Optical Modes in Room-Temperature Ionic Liquids.

The terahertz spectrum encodes information about dynamics, structure, and intermolecular interactions of liquids being probed both experimentally and computationally by techniques such as Raman and far-infrared spectroscopies and molecular dynamics (MD) simulation. In the case of room temperature ionic liquids (RTILs), there has been a debate whether a mode observed at about 1.5-2.7 THz (50-90 cm-1) is due to a quasi-lattice structure or the formation of complex ions. Here we show through the analysis of Raman and far-infrared spectra and MD simulation of a typical RTIL that this mode has a collective optic-like character. Then, employing a simple model based on the theory for optical phonons in crystals, we show that a correlation between the frequency of this mode and material parameters holds for different RTILs. These results, which encompass a wide range of samples, reinforce a quasi-lattice view of the liquid phase.

Ion pair free energy surface as a probe of ionic liquid structure.

Numerous combinations of cations and anions are possible for the production of ionic liquids with fine-tuned properties once the correlation with the molecular structure is known. In this sense, computer simulations are useful tools to explain and even predict the properties of ionic liquids. However, quantum mechanical methods are usually restricted to either small clusters or short time scales so that parameterized force fields are required to study the bulk liquids. In this work, a method is proposed to enable a comparison between the quantum mechanical system and both polarizable and nonpolarizable force fields by means of the calculation of free energy surfaces for the translational motion of the anion around the cation in gas phase. This method was tested for imidazolium-based cations with 3 different anions, [BF4]-, [N(CN)2]-, and [NTf2]-. Better agreement was found with the density functional theory calculations when polarizability is introduced in the force field. In addition, the ion pair free energy surfaces reproduced the main structural patterns observed in the first coordination shell in molecular dynamics simulations of the bulk liquid, proving to be useful probes for the liquid phase structure that can be computed with higher level methods and the comparison with forcefields can indicate further improvements in their parameterization.

Low-Temperature Phase Transitions of the Ionic Liquid 1-Ethyl-3-methylimidazolium Dicyanamide.

Several calorimetric measurements have shown that 1-ethyl-3-methylimidazolium dicyanamide, [C2C1im][N(CN)2], is a glass-forming liquid, even though it is a low-viscous liquid at room temperature. Here, we found slow crystallization during cooling of [C2C1im][N(CN)2] along Raman spectroscopy measurements. The low-frequency range of the Raman spectrum shows that the same crystalline phase is obtained at 210 K either by cooling or by reheating the glass (cold-crystallization). Another crystalline phase is formed at ca. 260 K just prior the melting at 270 K. X-ray diffraction and calorimetric measurements confirm that there are two crystalline phases of [C2C1im][N(CN)2]. The Raman spectra indicate that polymorphism is related to [C2C1im]+ with the ethyl chain on the plane of the imidazolium ring (the low-temperature crystal) or nonplanar (the high-temperature crystal). The structural reason for the glass-forming ability of [C2C1im][N(CN)2], despite the relatively simple molecular structures of the ions, was pursued by quantum chemistry calculations and molecular dynamics (MD) simulations. Density functional theory calculations were performed for ionic pairs in order to draw free-energy surfaces of the anion around the cation. The MD simulations using a polarizable model provided maps of occurrence of anions around cations. Both the quantum and classical calculations suggest that the delocalization of preferred positions of the anion around the cation, which adopts different conformations of the ethyl chain, is on the origin of the crystallization being hampered during cooling and the resulting glass-forming ability of [C2C1im][N(CN)2].

On the Regular Behavior of a Binary Mixture of Ionic Liquids.

Binary mixtures of the ionic liquids 1-ethyl-3-methylimidazolium dicyanamide and tributyl(methyl)ammonium dicyanamide, [C2C1im]x[N4441](1-x)[N(CN)2], are studied by means of their excess properties (enthalpy and volume), their viscosity, and their Raman spectra. The mixtures exhibit positive values of excess volume VE and excess enthalpy HE. The plot of the logarithm of viscosity as a function of composition is consistent with the finding HE > 0. The excess thermodynamic properties of the ionic liquid mixtures are compared with well-known results for binary mixtures of (high-temperature) molten salts. The asymmetry of the curve HE versus composition is reproduced by considering the volumetric fraction of the [C2C1im]x[N4441](1-x)[N(CN)2] mixtures according to basic results of the theory of regular solutions. The anion totally symmetric stretching mode νs(C≡N) was used to probe the local environment around the anion in a Raman spectroscopy investigation of the [C2C1im]x[N4441](1-x)[N(CN)2] mixtures as a function of the composition. The Raman band shape (peak frequency and bandwidth) of the νs(C≡N) mode changes as expected from the fluctuation of concentration around the probe oscillator according to models for vibrational spectroscopy of liquid mixtures. This thermodynamic and spectroscopic study leads to the classification of the [C2C1im]x[N4441](1-x)[N(CN)2] mixture as a regular solution. The application of the regular solution equation for HE to other ionic liquid binary mixtures is discussed.

Low-frequency Raman spectra of a glass-forming ionic liquid at low temperature and high pressure.

The frequency range below ∼100 cm-1 of the Raman spectrum of a glass-forming liquid exhibits two features that characterize the short-time (THz) dynamics: the quasi-elastic scattering (QES) tail and the boson peak (BP). In this work, we follow temperature and pressure effects on the intermolecular dynamics of a typical ionic liquid, 1-butyl-1-methylpiperidinium bis(trifluoromethanesulfonyl)imide, [Pip14][[NTf2]. The glass transition temperature of [Pip14][[NTf2] at atmospheric pressure is Tg = 198 K, and the pressure of glass transition at room temperature is Pg = 1.1 GPa. Raman spectra obtained while cooling the liquid or heating the glass exhibit hysteresis in QES and BP intensities, IQES and IBP. The dependence of IQES, IBP, and the BP frequency, ωBP, with pressure up to the glass transition is steeper than the temperature dependence due to the stronger pressure effect on density within the GPa range. The temperature and pressure behaviors of the parameters IQES, IBP, and ωBP obtained here for [Pip14][[NTf2] are discussed in light of known results for other glass-formers.

Spatial and thermal signatures of α and ß relaxations in glassy and glacial aliphatic ionic liquids.

The competition between Coulomb and van der Waals interactions brings forth unique dynamic features and broad applications to ionic liquids. Herein, we present a combined calorimetric, X-ray diffraction, incoherent elastic, and quasi-elastic neutron scattering study, over a wide temperature range (180-340 K), of the relaxational dynamics of the liquid, supercooled liquid, crystalline, glassy, and glacial states of two model ionic liquids: tributylmethylammonium (a good glass-former) and butyltrimethylammonium (a good crystal-former) cations and the bis(trifluoromethanesulfonyl)imide anion. In both systems, we observed two distinct relaxation processes. The Q-dependence of the respective relaxation time shows that the α-process is diffusive, while the ß-process is modulated by the structure of the liquids.

Structure and Reactivity of the Ionic Liquid 1-Allyl-3-methylimidazolium Iodide under High Pressure.

Poly(ionic liquid)s are an interesting class of compounds because of their unique chemical and physical properties gathering the characteristics of ionic liquids and polymers. Pressure and temperature have been demonstrated to be alternative parameters to obtain polymers from monomeric species using only physical tools. In this work, we investigate the reaction under high pressure and room temperature of the ionic liquid 1-allyl-3-methylimidazolium iodide by using the diamond anvil cell technique in combination with synchrotron X-ray diffraction and electronic and vibrational spectroscopies. The results indicate a chemical reaction happening through the terminal double bond of the allyl group both in crystalline and glassy phases with the onset of the reaction around ∼7 GPa. Vibrational spectra present evidence for an oligomerization reaction in both the phases. The reaction occurring both in glassy and crystal phases indicates a mechanism not driven by collective motions and likely related to local topological arrangements. The results presented herein extend our understanding of ionic liquid instability boundaries under high pressure and contribute to the development of alternative synthetic routes to achieve poly(ionic liquids).

Mechanical heterogeneity in ionic liquids.

Molecular dynamics (MD) simulations of five ionic liquids based on 1-alkyl-3-methylimidazolium cations, [C n C1im]+, have been performed in order to calculate high-frequency elastic moduli and to evaluate heterogeneity of local elastic moduli. The MD simulations of [C n C1im][NO3], n = 2, 4, 6, and 8, assessed the effect of domain segregation when the alkyl chain length increases, and [C8C1im][PF6] assessed the effect of strength of anion-cation interaction. Dispersion curves of excitation energies of longitudinal and transverse acoustic, LA and TA, modes were obtained from time correlation functions of mass currents at different wavevectors. High-frequency sound velocity of LA modes depends on the alkyl chain length, but sound velocity for TA modes does not. High-frequency bulk and shear moduli, K ∞ and G ∞ , depend on the alkyl chain length because of a density effect. Both K ∞ and G ∞ are strongly dependent on the anion. The calculation of local bulk and shear moduli was accomplished by performing bulk and shear deformations of the systems cooled to 0 K. The simulations showed a clear connection between structural and elastic modulus heterogeneities. The development of nano-heterogeneous structure with increasing length of the alkyl chain in [C n C1im][NO3] implies lower values for local bulk and shear moduli in the non-polar domains. The mean value and the standard deviations of distributions of local elastic moduli decrease when [NO3]- is replaced by the less coordinating [PF6]- anion.

Communication: Glass transition and melting lines of an ionic liquid.

The phase diagram of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethanesufonyl)imide, [Pyrr1,4][NTf2], was explored by synchroton X-ray diffraction and Raman scattering measurements as a function of temperature and pressure. Glass transition Tg(p) and melting Tm(p) temperatures were obtained from atmospheric pressure up to ca. 2.0 GPa. We found that both the Tg(p) and Tm(p) curves follow essentially the same pressure dependence. The similarity of pressure coefficients, dTg/dp ≈ dTm/dp, is explained within the non-equilibrium thermodynamics approach for the glass transition by assuming that one of the Ehrenfest equations is appropriated for Tg(p), whereas Tm(p) follows the Clausius-Clapeyron equation valid for the first-order transitions. The results highlight that ionic liquids are excellent model systems to address fundamental questions related to the glass transition.

Ultraslow Phase Transitions in an Anion-Anion Hydrogen-Bonded Ionic Liquid.

A Raman spectroscopy study of 1-ethyl-3-methylimidazolium hydrogen sulfate, [C2C1im][HSO4], as a function of temperature, has been performed to reveal the role played by anion-anion hydrogen bond on the phase transitions of this ionic liquid. Anion-anion hydrogen bonding implies high viscosity, good glass-forming ability, and also moderate fragility of [C2C1im][HSO4] in comparison with other ionic liquids. Heating [C2C1im][HSO4] from the glassy phase results in cold crystallization at ∼245 K. A solid-solid transition (crystal I → crystal II) is barely discernible in calorimetric measurements at typical heating rates, but it is clearly revealed by Raman spectroscopy and X-ray diffraction. Raman spectroscopy indicates that crystal I has extended ([HSO4]-)n chains of hydrogen-bonded anions but crystal II has not. Raman spectra recorded at isothermal condition show the ultraslow dynamics of cold crystallization, solid-solid transition, and continuous melting of [C2C1im][HSO4]. A brief comparison is also provided between [C2C1im][HSO4] and [C4C1im][HSO4], as Raman spectroscopy shows that the latter does not form the crystalline phase with extended anion-anion chains.

Structure and Liquid Fragility in Sodium Carbonate.

The relationship between local structure and dynamics is explored for molten sodium carbonate. A flexible fluctuating-charge model, which allows for changes in the shape and charge distribution of the carbonate molecular anion, is developed. The system shows the evolution of highly temperature-dependent complex low-dimensional structures which control the dynamics (and hence the liquid fragility). By varying the molecular anion charge distribution, the key interactions responsible for the formation of these structures can be identified and rationalized. An increase in the mean charge separation within the carbonate ions increases the connectivity of the emerging structures and leads to an increase in the system fragility.